The  Astrophysical  Journal  Supplement  Series,  155:89-100,  2004  November 

©  2004.  The  American  Astronomical  Society.  All  rights  reserved.  Printed  in  U.S.A. 


A  wide-field,  low-frequency  radio  survey  of  the  field  of  M3 1.  I. 

CONSTRUCTION  AND  STATISTICAL  ANALYSIS  OF  THE  SOURCE  CATALOG 

Joseph  D.  Gelfand 

Harvard-Smithsonian  Center  for  Astrophysics,  MS- 10,  60  Garden  Street,  Cambridge,  MA  02138;  jgelfand@cfa.harvard.edu 

T.  Joseph  W.  Lazio 

Naval  Research  Laboratory,  Code  7213,  4555  Overlook  Avenue  SW,  Washington,  DC  20375-5351;  joseph.lazio@nrl.navy.mil 

AND 

B.  M.  Gaensler 

Harvard-Smithsonian  Center  for  Astrophysics,  60  Garden  Street,  Cambridge,  MA  02138;  bgaensler@cfa.harvard.edu 
Received  2004  January  24;  accepted  2004  April  27 


ABSTRACT 

We  present  here  the  results  of  a  325  MHz  radio  survey  of  M31,  conducted  with  the  A  configuration  of  the  Very 
Large  Array.  The  survey  covered  an  area  of  7.6  deg^,  and  a  total  of  405  radio  sources  between  ^6"  and  170"  in 
extent  were  mapped  with  a  resolution  of  6"  and  a  1  cr  sensitivity  of  ^0.6  mJy  beam^F  For  each  source,  its 
morphological  class,  major  axis  9m,  minor  axis  9„,  position  angle  0pA,  peak  flux  7,  integrated  flux  density  S, 
spectral  index  a,  and  spectral  curvature  parameter  ^p  were  calculated.  A  comparison  of  the  flux  and  radial 
distribution — both  in  the  plane  of  the  sky  and  in  the  plane  of  M31 — of  these  sources  with  those  of  the  XMM- 
Newton  Large-Scale  Structure  Survey  and  Westerbork  Northern  Sky  Survey  revealed  that  a  vast  majority  of 
sources  detected  are  background  radio  galaxies.  As  a  result  of  this  analysis,  we  expect  that  only  a  few  sources  are 
intrinsic  to  M3 1 .  These  sources  are  identified  and  discussed  in  an  accompanying  paper. 

Subject  headings:  catalogs  —  galaxies:  individual  (M31)  —  radio  continuum:  galaxies  — 
radio  continuum:  general 

Online  material:  additional  figures,  machine-readable  table 


1.  INTRODUCTION 

In  order  to  understand  the  dynamics  of  a  galaxy  it  is  crucial 
to  know  properties  of  its  radio  population.  Much  has  been 
learned  about  the  Milky  Way  from  radio  observations  of  the 
supernova  remnants  (SNRs),  H  ii  regions,  and  pulsars  detected 
within  it.  However,  a  complete  census  of  radio  sources  in  the 
Milky  Way  is  difficult  to  obtain  because  of  source  confusion. 
In  addition,  since  the  distances  to  these  sources  are  often 
extremely  uncertain,  it  is  difficult  to  determine  their  proper¬ 
ties.  As  a  result,  it  is  advantageous  to  observe  external  gal¬ 
axies  to  learn  about  the  dynamics  and  properties  of  “normal” 
galaxies.  The  first  step  of  this  process  is  to  obtain  a  census  of 
radio  populations  in  a  galaxy,  and  it  is  for  this  reason  that  we 
have  surveyed  the  radio  population  of  M3 1 ,  the  nearest  spiral 
galaxy. 

M3 1  has  been  surveyed  in  the  radio  before,  both  as  parts  of 
larger  surveys  such  as  the  Westerbork  Northern  Sky  Survey 
( WENSS;  Rengelink  et  al.  1997)  and  NRAO  VLA  Sky  Survey 
(NVSS;  Condon  et  al.  1998),  and  as  the  focus  of  dedicated 
surveys  such  as  the  36W  (Bystedt  et  al.  1984),  37W  (Walterbos 
et  al.  1985),  and  Braun  (1990)  surveys,  the  properties  of  all  of 
which  are  summarized  in  Table  1 .  Despite  its  proximity  to  the 
Milky  Way,  mapping  M3 1  is  difficult  because  its  large  angular 
size  on  the  sky  (>2°)  requires  many  pointings  at  higher  radio 
frequencies  GHz)  to  fully  cover.  As  a  result,  existing 

surveys  of  M3 1  are  either  deep  but  cover  only  a  small  region 
of  M3 1 ,  or  cover  the  entire  optical  disk  of  M3 1  with  relatively 
poor  sensitivity,  as  shown  in  Table  1  and  Figure  1.  Because  of 
this,  the  radio  population  of  discrete  sources  in  M3 1  is  not  well 
understood.  To  rectify  this  situation,  we  have  surveyed  M31 


with  the  Very  Large  Array'  (VLA)  at  325  MHz  (1  =  90  cm) 
using  the  A  configuration,  achieving  a  1  cr  sensitivity  less  than 
a  mJy  with  a  resolution  of  6"  (^20  pc  at  the  distance  of  M3 1 , 
assumed  to  be  780  kpc;  Stanek  &  Gamavich  1998)  over  the 
entire  optical  disk  of  M31  because  of  the  large  size  of  the 
VLA’s  primary  beam  at  this  frequency.  Using  the  A  configu¬ 
ration  has  the  advantage  that  extended  emission  from  M3 1  was 
resolved  out,  allowing  us  to  better  determine  the  properties 
of  the  compact  radio  population.  As  seen  in  Table  1,  the  survey 
presented  here  does  very  well  in  resolution  {9ies)  versus  field 
of  view  (FOV)  when  compared  to  past  surveys  of  M31  and 
has  a  similar  sensitivity  to  that  of  previous  higher  frequency 
surveys. 

This  paper  describes  the  observations,  the  data  reduction 
process  that  led  to  the  final  source  list,  and  the  statistical 
properties  of  the  detected  sources.  This  paper  is  structured 
as  follows:  §  2  describes  the  observations  and  data  reduction 
techniques  used  to  make  the  final  source  list,  §  3  presents  the 
final  (GLG)  source  list,  §  4  presents  the  statistical  properties 
of  the  GLG  source  list,  and  §  5  presents  the  conclusions 
derived  from  this  analysis.  In  a  companion  paper  (J.  D. 
Gelfand,  T.  J.  W.  Lazio,  &  B.  M.  Gaensler  2004,  in  prep¬ 
aration),  we  will  classify  the  sources  and  discuss  their 
properties. 

This  paper  uses  the  convention  of  is  the  flux  density  of  a 
source  at  a  frequency  u  MHz,  radio  spectral  index  a  is  defined 

*  The  VLA  is  operated  by  The  National  Radio  Astronomy  Observatory, 
which  is  a  facility  of  the  National  Science  Foundation  operated  under  coop¬ 
erative  agreement  by  Associated  Universities,  Inc. 


89 


Report  Documentation  Page 

Form  Approved 

0MB  No.  0704-0188 

Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 

VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  0MB  control  number. 

1.  REPORT  DATE 

NOV  2004 

2.  REPORT  TYPE 

3.  DATES  COVERED 

00-00-2004  to  00-00-2004 

4.  TITLE  AND  SUBTITLE 

5a.  CONTRACT  NUMBER 

A  Wide-Field,  Low-Frequency  Radio  Survey  of  the  Field  of  M31. 1. 
Construction  and  Statistical  Analysis  of  the  Source  Catalog 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Naval  Research  Laboratory, Code  7213,4555  Overlook  Avenue, 

SW, Washington, DC, 20375 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

10.  SPONSOR/MONITOR’S  ACRONYM(S) 

11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

15.  SUBJECT  TERMS 

16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 
ABSTRACT 

18.  NUMBER 

OF  PAGES 

12 

19a.  NAME  OF 
RESPONSIBLE  PERSON 

a.  REPORT 

unclassified 

b.  ABSTRACT 

unclassified 

c.  THIS  PAGE 

unclassified 

Standard  Form  298  (Rev.  8-98} 

Prescribed  by  ANSI  Std  Z39-18 


90 


GELFAND,  LAZIO,  &  GAENSLER 


Vol.  155 


TABLE  1 

Summary  of  Radio  Surveys 


Flux  Density 

Observation  Frequency  Limit^  Resolution  FOV  Number  of 


Survey  Date  (MHz)  (niJy)  (arcsec)  (deg^)  Sources  Reference 


WENSS .  1991-1993  325  18  54  x  82  ...  288  Rengelink  et  al.  (1997) 

36W .  1975  610  3  54  x  82  10.1  392  Bystedt  et  al.  (1984) 

37W .  1975-1979  1412  1  23  x  35  3.1  249  Walterbos  et  al.  (1985) 

Braun .  1990  1465  0.15  5  1.9  534  Braun  (1990) 

NVSS .  1993-1996  1400  2.5  45  ...  831  Condon  et  al.  (1998) 

GLG .  2000  325  3  6  7.6  405  This  paper 


^  The  flux  density  limits  given  for  WENSS,  36W,  37W,  Braun,  and  NVSS  surveys  are  the  lowest  flux  densities  quoted  by  these  authors.  For  all 
surveys  but  the  36W  survey,  this  corresponds  to  a  5  (T  cutoff.  For  the  36W  survey,  this  corresponds  to  a  4  cr  cutoff. 


as  Si,  fx.v°‘ ,  and  the  distance  to  M3 1  is  assumed  to  be 
780  ±  13  kpc  (Stanek  &  Gamavich  1998). 

2.  OBSERVATIONS  AND  DATA  REDUCTION 

This  section  describes  the  observations  and  data  reduction 
process  that  led  to  the  GLG  source  list.  As  detailed  below,  the 
data  reduction  process  dilFers  from  that  of  higher  frequency 
observations  because  of  the  large  field  of  view.  In  addition, 
325  MHz  observations  are  plagued  by  radio  frequency  inter¬ 
ference  (RFI),  which  must  be  removed  before  imaging. 

2.1.  Observations 

This  paper  is  the  result  of  a  5  hr  {~4  hr  on-source)  obser¬ 
vation  of  M3 1  conducted  on  2000  December  1 5  with  the  VLA 
A  configuration.  The  observations  were  done  in  spectral  line 
mode  with  two  intermediate  frequencies  (IFs),  centered  at 
321.6  and  327.5  MHz,  each  with  a  bandwidth  of  3125  kHz 
divided  into  64  frequency  channels.  Spectral  line  mode  was 
used  to  reduce  the  effect  of  primary  beam  chromatic  aberration 
(bandwidth  smearing)  and,  more  importantly,  to  expedite  ex¬ 
cision  of  radio  frequency  interference  (RFI).  The  receivers 
recorded  both  right  circular  polarization  (RCP)  and  left  circular 
polarization  (LCP)  data,  which  means  that  only  Stokes  I  (total 
intensity)  and  Stokes  V  (circular  polarization)  were  measured. 
Calibrator  sources  3C  405  (Cyg  A)  and  TXS  0035+413  were 
observed  hourly  for  ~5  minutes.  Cyg  A  was  used  to  set  the  flux 
density  scale  and  calibrate  the  bandpass  by  applying  a  333  MHz 
model  to  the  observed  image. ^  This  model  assumes  that  Cyg  A 
has  a  total  flux  density  of  5983  Jy  at  325  MHz,  within  the 
6%  error  of  the  Baars  flux  density  scale  (Baars  et  al.  1977). 
The  calibrator  source  TX  0035+413  was  then  used  to  cali¬ 
brate  the  visibilities.  Table  2  summarizes  the  properties  of  the 
observation.^ 

2.2.  Data  Reduction 

The  initial  phase  and  flux  density  calibration  was  done  using 
the  AIPS  task  CALIB.  For  both  Cyg  A  and  TXS  0035+413, 
edge  channels  1-5  and  58-63  were  flagged  owing  to  poor  fre¬ 
quency  response  and  the  visibility  data  were  inspected  manu¬ 
ally  for  RFI.  During  the  first  half  of  the  first  scan  of  Cyg  A, 
antenna  1 3  malfunctioned.  To  preserve  data  from  this  antenna 

^  Available  at  http://rsd-www.nrl. navy.mil/7213/lazio/tutorialA^LAmodels/ 
models.html. 

^  The  FOV  listed  in  this  table  is  equal  to  FOV  —  where  R  is  the 
distance  in  degrees  from  the  pointing  center  to  the  most  distant  source  detected. 
This  is  larger  than  the  FWHM  of  the  primary  beam  owing  to  the  shapes  of  the 
“facets”  discussed  in  §  2.2.1. 


taken  in  the  second  half  of  this  scan,  all  data  points  involving 
antenna  13  in  the  first  half  of  this  scan  were  flagged  and  the 
phase  and  flux  solutions  calculated  for  the  second  half  were 
renormalized.  After  the  flux  and  phase  calibration  solutions 
were  calculated,  they  were  applied  to  all  the  visibility  data, 
including  the  M31  visibilities. 

The  process  of  converting  the  raw  data  into  a  usable  image 
and  source  list  was  a  long  and  complicated  one,  as  illustrated 
in  Figure  2.  The  major  complication  in  this  process  is  the 
presence  of  RFI,  and  it  was  necessary  to  remove  as  much  of  it 
as  possible  before  any  imaging.  This  was  accomplished  by 
inspecting  each  baseline  visually  for  RFI.  Both  the  RCP  and 
LCP  data  were  searched  separately  for  RFI  because,  being 
terrestrial  in  origin,  it  was  often  polarized  and  would  appear  in 
one  polarization  channel  but  not  the  other.  To  avoid  compli¬ 
cations  when  calculating  the  total  intensity  (Stokes  I),  flags  in 
one  polarization  were  applied  to  the  other.  Since  polarized  RFI 
is  strong  in  Stokes  V  (Stokes  V  =  RCP— LCP),  all  baselines 
were  searched  for  RFI  in  this  Stokes  parameter  as  well. 

2.2.1.  Image  Production 

Only  after  all  the  RFI  was  removed  was  it  possible  to  pro¬ 
duce  a  high-quality  image  of  M3 1 .  Owing  to  its  large  FOV  at 
325  MHz,  the  VLA  is  not  coplanar  at  this  frequency.  Therefore, 
in  order  to  make  the  image  using  standard  radio  techniques,  the 
primary  beam  was  divided  into  73  small  “facets”  over  which 
the  VLA  can  be  safely  assumed  to  be  coplanar  (Cornwell  & 
Perley  1992). 

Two  sets  of  images  were  made,  both  created  by  an  iterative 
process  of  CLEANing  and  self-calibration.  For  Images  A,  the 
self-calibration  process  used  every  positive  CLEAN  compo¬ 
nent  (CC).  As  a  result,  the  calibration  model  included  CCs  that 
were  peaks  in  the  noise,  artificially  lowering  the  noise  of  the 
image  and  increasing  the  probability  of  a  false  detection. 
However,  these  maps  were  good  enough  to  find  sources 
(Source  list  I)  that  generated  the  restrictive  clean  boxes 
(CLEAN  Boxes  2)  used  to  image  Images  B. 

The  second  set  of  images  (Images  B),  which  generated  the 
source  list  presented  in  this  paper,  were  also  created  using  an 
iterative  process  involving  the  CLEAN  algorithm  and  self¬ 
calibration.  However,  this  time  CCs  were  only  placed  around 
sources  in  Source  list  I  and  off-axis  WENSS  sources  (CLEAN 
Boxes  2).  The  use  of  smaller  clean  boxes  decreased  the  pos¬ 
sibility  of  “overcleaning”  and  removes  any  CLEAN  bias  (see 
Condon  et  al.  1998  for  a  good  description  of  CLEAN  bias).  In 
addition,  the  self-calibration  process  used  only  the  strongest 
CCs  inside  the  primary  beam  to  maximize  the  accuracy  of  the 
phase  calibration  model  inside  the  primary  beam.  Using  this 
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Fig.  1. — Contours  showing  the  1  cr  sensitivity  of  the  (clockwise  from  upper  left)  GLG  (this  paper),  36W  (Bystedt  et  al.  1984),  Braun  (1990),  and  37W  (Walterbos 
et  al.  1985)  radio  surveys  across  the  field  of  M31.  The  sensitivity  for  the  36W,  Braun,  and  37W  surveys  were  generated  using  the  procedure  described  in  §  4.3.  The 
optical  image  is  of  M3 1  from  the  Palomar  Optical  Sky  Survey.  For  the  GLG  and  36W  sensitivity  maps,  the  contour  levels  are  0.6,  0.65,  0.7,  . . . ,  1.0,  1.25, 1.5 , ... , 
3.0  mjy  beam“’.  The  jagged  edges  of  the  GLG  sensitivity  map  are  due  to  the  shapes  of  the  facets  discussed  in  §  2.2.1.  For  the  37W  sensitivity  map,  the  contour 
levels  are  0.1,  0.15,  0.2,  ...,  1.0  mJy  beam"',  and  for  the  Braun  sensitivity  map,  the  contour  levels  are  0.025,  0.035,  0.045,  0.05,  0.06, . . . ,  0.1,  0.25  mJy  beam"*. 


method,  we  generated  maps  with  a  typical  noise  level  of 
~0.6  mJy  beam^*  and  a  signal-to-noise  ratio  of  ~1000  around 
the  strongest  sources — equivalent  to  that  of  the  highest  quality 
low-frequency  radio  images. 

2.2.2.  Generating  the  Source  List 

Sources  in  the  first  (Images  A)  and  second  (Images  B)  set  of 
images  were  found  using  the  MIRIAD  task  SFIND  (ver.  2.0), 
which  uses  a  false  discovery  rate  (FDR)  algorithm  to  detect 
sources  (Hopkins  et  al.  2002).  This  algorithm  is  similar  to  a 
Gaussian  source-finding  algorithm,  except  that  it  also  accounts 
for  noise  variations  on  a  user-defined  size  (the  rmsbox  pa¬ 
rameter).  For  each  one  of  the  central  73  facets,  we  made  six 


source  lists.  The  first  three  source  lists  were  made  using 
rmsbox  =  1.25,  1.88,  and  2.5  beams  (20,  30,  and  40  pixels). 
Multiple  values  were  used  for  the  rmsbox  parameter  because 
sources  the  size  of  rmsbox  were  often  missed  by  SFIND,  and 
values  of  rmsbox  larger  than  2.5  beams  caused  SFIND  to  miss 
many  sources.  For  the  different  values  of  rmsbox,  SFIND 
detected  a  similar  number  of  sources  but  there  were  sources 
found  for  one  value  but  not  others.  These  three  source  lists 
were  then  compiled  into  one  list.  Since  the  facets  overlapped 
slightly,  sources  on  the  edge  of  a  facet  were  “discovered”  more 
than  once.  When  merging  the  source  list  of  each  facet  into  a 
single  source  list,  the  source  with  the  strongest  signal-to-noise 
from  this  set  of  duplicates  was  kept.  SFIND  also  had  difficulty 
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TABLE  2 

Observation  Properties 

A.  Observed  Sources 

Integration  Time 

Source 

(minutes) 

Number  of  Visibilities 

3C  405  . 

25.8 

47330 

TXS  0035+413  . 

23.1 

42250 

M3 1 . 

239.5 

438657 

B.  Receiver  Properties 

Frequency 

Bandwidth  (Channel  Separation) 

IF  Number 

(MHz) 

(kHz) 

1 . 

327.50 

3125.0  (48.83) 

2 . 

321.56 

3125.0  (48.83) 

C.  Primary  Beam  Properties 

Property 

Value 

Full  width  at  half-maximum 

138' 

Field  of  view . 

7.6  deg^ 

Resolution . 

6"  X  6" 

rms  at  center  of  field . 

0.61  mjy 

D.  General  Properties 

Property 

Value 

Date  of  observation . 

2000  Dec  15 

Pointing  center  (J2000) . 

R.A.  =  00‘'42”’45t9,  deck  =  41°16T8"0 

with  large  and/or  complicated  sources,  the  integrated  proper¬ 
ties  of  these  sources  were  determined  by  running  SFIND 
without  its  FDR  capabilities"^  on  maps  smoothed  to  a  resolution 
of  20"  X  20".  In  these  smoothed  maps,  a  complicated  source 
(e.g.,  a  resolved  radio  jet)  would  be  unresolved,  allowing 
SFIND  to  better  determine  its  integrated  properties.  As  with  the 
high-resolution  maps,  three  source  lists  were  made,  each  with 
a  different  value  of  rmsbox,  rmsbox  =  3.13,  6.25,  and  9.38 
beams  (50,  100,  and  150  pixels).  Again,  these  three  source  lists 
were  compiled  into  one.  Once  the  final  high-resolution  and 
smoother  source  lists  were  made,  they  were  compiled  into  one. 
In  the  case  where  a  source  in  the  smoothed  map  overlapped 
many  sources  in  the  high-resolution  maps,  the  smoothed 
source  was  kept.  Sources  that  appeared  only  in  the  smoothed 
maps  were  kept  as  well.  However,  in  all  other  cases  the  source 
in  the  high-resolution  maps  were  kept.  The  positions  and  sizes 
of  sources  found  in  Images  A  (Source  list  I)  defined  the  clean 
boxes  (CLEAN  Boxes  2)  used  to  make  the  second  and  final  set 
of  images  (Images  B).  The  final  source  list  presented  here 
(hereafter  the  GLG  source  list),  given  in  Table  3,  was  made 
from  Images  B  using  this  procedure.  The  location  of  these 
sources  is  shown  in  Figure  3. 

2.3.  Source  Properties 

The  MIRIAD  task  SFIND  determines  the  R.A.,  deck,  local 
flux  density  rms  a^s,  peak  brightness  I,  integrated  flux  density 
S,  image  major  axis  9m,  image  minor  axis  9„,  and  position 
angle  0pA  of  every  source  it  finds.  Sources  were  separated  into 

Without  the  FDR  algorithm,  SFIND  works  like  other  Gaussian  source 
finding  programs — looking  for  regions  with  a  signal-to-noise  ratio  greater  than 
a  user-defined  value.  In  our  case,  we  instructed  SFIND  to  keep  all  >5  u 
sources. 


four  morphological  categories — unresolved  (U),  elongated 
(El),  complex  (C),  and  extended  (Ex) — defined  in  the  fol¬ 
lowing  way: 

1.  Unresolved  sources  are  sources  in  the  high-resolution 
image  with  9m  <2x9^. 

2.  Elongated  sources  are  also  sources  in  the  high-resolution 
image,  but  with  0m  >  2  x  0„,.  In  practice,  the  major  distinction 
between  unresolved  and  elongated  sources  is  that  the  major 
axis  is  resolved  for  elongated  sources. 

3.  Complex  sources  are  sources  in  the  smoothed  maps  that 
overlap  several  sources  in  the  high-resolution  maps. 

4.  Extended  sources  are  sources  in  the  smoothed  maps  that 
do  not  overlap  sources  in  the  high-resolution  maps.  We  defined 
“ExG”  sources  to  be  extended  sources  detected  in  other  radio 
surveys  of  M3 1 .  Extended  sources  without  counterparts  in  other 
radio  surveys  (“Ex”  sources)  are  most  likely  false  detections  but 
are  included  here  for  completeness. 

A  total  of  405  sources  were  detected:  281  are  unresolved,  16 
elongated,  51  complex,  and  57  extended — of  which  two  are 
ExG  sources. 

The  source  flux  density  S  was  calculated  by  correcting  the 
value  returned  by  SFIND  for  primary  beam  attenuation.  After 
the  primary  beam  correction  was  applied,  the  errors  in  9m,  9„„ 
and  0pA  of  each  source  were  calculated  using  equations  (29a), 
(29b),  and  (30)  in  Condon  et  al.  (1998).  In  this  calculation,  we 
assumed  that  the  signal-to-noise  of  a  source  is  equal  to  S/tXfnjs. 
We  also  calculated  the  error  in  integrated  flux  density  (S) 
and  peak  brightness  (/)  using  equations  (36a)  and  (37)  in 
Condon  et  al.  (1998)  for  a  resolved  (C,  Ex,  and  ExG)  source, 
equations  (39)  and  (41)  in  Condon  et  al.  (1998)  for  a  source 
where  only  the  major  axis  is  resolved  (El),  and  equation  (43)  in 
Condon  et  al.  (1998)  for  an  unresolved  (U)  source.  However, 
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TABLE  3 
GLG  Source  List 


Name 

Type 

R.A. 

(arcsec) 

Decl. 

0-6 

(arcsec) 

^M31 

(arcmin) 

^M31 

(kpc) 

Sm 

(arcsec) 

<y(8M) 

Om 

(arcsec) 

^PA 

(deg) 

<^PA 

/ 

(mJy  beam"*) 

OA 

5 

(mJy) 

<rs 

*^nns 

a 

GLG  001 . 

U 

00  42  50.40 

+41  16  10.28 

0.8 

5.7 

4.6 

15 

7.94 

8.26 

0.69 

-0.64  ±  0.06 

0.34 

0.40 

0.7 

0.4 

0.3 

6 

0.70 

0.70 

GLG  002 . 

U 

00  42  26.62 

+41  18  04.62 

4.0 

4.9 

4.5 

76 

4.17 

3.65 

0.63 

-1.45  ±  0.34 

0.71 

0.66 

3.9 

1.4 

1.2 

155 

0.59 

0.59 

0.13  ±  0.05 

GLG  003 . 

u 

00  42  32.90 

+41  11  43.57 

5.1 

5.0 

4.7 

14 

5.12 

4.76 

0.54 

-0.95  ±  0.02 

0.46 

0.48 

1.4 

0.6 

0.5 

26 

0.53 

0.53 

GLG  004 . 

u 

00  43  18.00 

+41  17  21.94 

6.1 

10.0 

7.7 

67 

131 

22.85 

0.60 

-0.84  ±  0.14 

0.22 

0.19 

4.2 

0.1 

0.0 

6 

0.70 

0.70 

0.04  ±  0.02 

GLG  005 . 

u 

00  43  26.34 

+41  19  12.04 

8.1 

5.7 

5.3 

12 

6.60 

7.97 

0.67 

<-1.15 

0.38 

0.41 

4.5 

0.4 

0.4 

14 

0.62 

0.62 

Notes. — Units  of  right  ascension  are  hours,  minutes,  and  seconds,  and  units  of  declination  are  degrees,  arcminutes,  and  arcseconds.  Table  3  is  available  in  its 
entirety  in  the  electronic  edition  of  the  Astrophysical  Journal  Supplement.  A  portion  is  shown  here  for  guidance  regarding  its  form  and  content. 


in  all  cases  we  set  the  “clean  bias”  term  {gb  in  the  NVSS 
paper;  Condon  et  al.  1998)  to  0  because  restrictive  clean  boxes 
were  used. 

To  check  the  flux  calibration,  we  compared  the  flux  densities 
of  GLG  sources  with  the  flux  densities  of  their  counterparts  in 
the  WENSS  catalog  (Rengelink  et  al.  1997),  which  was  also 
conducted  at  325  MHz.  Figure  4  compares  the  WENSS  and 
GLG  flux  density  for  all  GLG  sources  with  a  “good”  match  in 
the  WENSS  catalog,  determined  using  the  criteria  discussed 
below  in  §  4.3.  In  general,  there  is  good  agreement  between  the 
flux  densities  measured  by  the  WENSS  and  GLG  surveys, 
particularly  for  the  strongest  sources  and  unresolved  sources — 
implying  that  the  flux  calibration  is  correct. 

After  the  flux  errors  were  calculated,  the  positions  of  the 
GLG  sources  were  registered  to  those  of  the  strongest  NVSS 


45“  40“  35“ 

Right  Ascension  (J2000) 


Fig.  3. — Palomar  Optical  Sky  Survey  image  of  M31  overlaid  with  the 
position  of  the  GLG  sources  {white  contours,  representing  the  major  and 
minor  axis  of  the  sources  given  in  Table  3)  and  the  1  a  sensitivity  of  the  GLG 
survey  {black  contours).  The  sensitivity  contours  correspond  to  0.6,  0.65,  0.7, 
0.75,  0.8,  0.85,  0.9,  0.95,  1.0,  1.5,  2.0,  and  2.5  mJy  beam-'. 


sources.  This  was  done  by  first  finding  all  GLG  sources  that 
overlapped  a  NVSS  source  with  5i4oo  >  15  mJy,  because  the 
positions  of  these  NVSS  sources  is  known  to  better  than  1" 
and  match  those  in  the  optical  reference  frame  to  ^40  mas 
(Condon  et  al.  1998).  Figure  5  shows  the  olFsets  between  the 
positions  of  these  15  GLG  sources  and  their  NVSS  counter¬ 
parts.  The  average  R.A.  offset  between  the  GLG  and  NVSS 
sources  was  — 0"099  ±  0^067,  and  the  average  deck  oflFset  was 
0"115±0"067,  the  error  being  the  standard  error  in  the 
mean.  Every  source  in  the  GLG  source  list  was  then  shifted  by 
the  average  offset,  and  the  intrinsic  uncertainty  in  R.A.  (cq.) 
and  deck  (e^)  were  set  to  the  error  of  the  offsets.  Errors  in 
the  positions  of  the  GLG  sources  were  then  calculated  with 
the  formulae  used  for  the  NVSS  source  list  (Condon  et  al. 
1998). 

3.  SOURCE  LIST 

The  final  source  list  is  presented  in  Table  3.  The  columns 
are  as  follows: 

Column  (I). — Name  of  the  source.  The  sources  are  ordered 
in  increasing  distance  from  the  center  of  the  FOV. 


10  100  1000  10000 
GLG  flux  [mJy] 


Fig.  4. — Flux  density  of  WENSS  sources  with  a  “good”  match  (as  defined 
in  §  4.3)  in  the  GLG  catalog  vs.  flux  density  of  corresponding  GLG  sources. 
The  diamonds  indicate  an  “U”  or  “El”  GLG  source,  and  the  sloping  straight 
line  indicates  where  the  WENSS  flux  density  equals  the  GLG  flux  density. 
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Fig.  5. — Positional  offset  (GLG-NVSS)  between  the  15  NVSS  and  GLG 
sources  with  S1400  >15  mJy.  These  sources  were  used  to  register  the  location 
of  the  GLG  sources. 

Column  (2). — Type  of  source.  This  is  the  morphological 
category  of  the  sources  as  defined  in  §  2.3.  The  abbreviation 
“U”  stands  for  unresolved,  “El”  for  elongated,  “C”  for 
complex,  “Ex”  for  extended,  and  “ExG”  for  an  extended 
source  with  a  counterpart  in  one  or  more  other  radio  surveys 
of  M31. 

Column  (3). — Right  ascension  of  the  source  (J2000).  Error 
in  R.A.  (cr„)  in  arcseconds. 

Column  (4). — Declination  of  the  source  (J2000).  Error  in 
decl.  {as)  in  arcseconds. 

Column  (5). — 0m31>  angular  distance  from  the  center  of 
M31,  in  arcmin.  Rm31>  projected  radius  of  source  in  M31  as 
defined  in  §  4.1,  in  kpc. 

Column  (6). — Major  axis  {9m)  of  the  source,  in  arcseconds. 
Error  in  major  axis  [cr(0M)]  in  arcseconds. 

Column  (7). — Minor  axis  {9^)  of  the  source,  in  arcseconds. 
Error  in  minor  axis  [a{9m)]  in  arcseconds. 

Column  (8). — Position  angle  (0pa)  north  through  east  of  the 
source,  in  degrees.  Error  in  position  angle  (ctpa)  in  degrees. 

Column  (9). — Peak  brightness  (7)  of  the  source,  in  mJy 
beam^*.  Error  in  peak  brightness  (ct/)  in  mJy  beam^*. 

Column  (10). — Integrated  flux  density  {S)  of  the  source  in 
mJy.  Error  in  integrated  flux  density  (<75),  in  mJy. 

Column  (11). — Root  mean  square  of  the  flux  density  around 
the  source  (tTrms),  in  mJy. 

Column  (12). — Radio  spectral  index  (a)  of  source,  as  cal¬ 
culated  in  §  4.3.  Spectral  curvature  parameter  (1^),  as  calculated 
in  §  4.3.  The  spectral  indices  of  Ex  sources  without  a  counter¬ 
part  in  another  radio  survey  were  not  published  because  this 
source  is  most  likely  spurious,  as  mentioned  in  §  2.3.  It  is  also 
possible  that  some  of  the  steeper  spectrum  U  sources  {a  <  —2) 
are  spurious  as  well. 

4.  PROPERTIES  OF  OBSERVED  SOURCES 


those  elsewhere,  and  if  so,  how  they  differed.  We  analyzed 
both  the  radial,  flux  density,  and  spectral  distribution  of  the 
GLG  sources,  and  for  comparison,  data  from  the  WENSS 
(Rengelink  et  al.  1997)  and  XMM-Newton  Large-Scale  Struc¬ 
ture  (VMM-LSS)  radio  survey  (Cohen  et  al.  2003)  surveys.  The 
WENSS  survey  mapped  the  entire  sky  north  of  b  =  30°  at 
V  =  325  MHz  with  a  limiting  flux  density  of  1 8  mJy  and  res¬ 
olution  of  54"  X  82"  around  M31  (Rengelink  et  al.  1997), 
while  the  325  MHz  field  of  the  XMM-LSS  radio  survey 
mapped  a  ^5.6  deg^  region  of  sky  with  a  resolution  of  6"3  and 
a  limiting  fiux  density  of  4  mJy  beam^^ — properties  similar  to 
the  observation  presented  here  (Cohen  et  al.  2003). 

4. 1 .  Radial  Distribution  of  Observed  Sources 

If  a  substantial  number  of  sources  intrinsic  to  M31  were 
detected,  we  would  expect  to  see  an  overdensity  of  sources  in 
the  optical  disk  of  M3 1 .  In  addition,  any  substructure  present 
in  M31  is  most  likely  symmetric  to  some  degree  and  might 
stand  out  in  the  radial  distribution  of  sources.  Unfortunately, 
interpreting  the  raw  radial  distribution  of  sources  is  dilficult 
because  the  gain  of  the  telescope  decreases  toward  the  edge  of 
the  FOV,  as  illustrated  in  Figure  1.  As  a  result,  faint  sources 
detectable  in  the  center  of  the  FOV  are  undetectable  at  the 
edges,  complicating  models  for  the  expected  observed  back¬ 
ground  radial  distribution  of  sources.  However,  since  the  XMM- 
LSS)  325  MHz  observation  was  also  done  using  one  pointing 
of  the  VLA  A  configuration,  its  primary  beam  shape  is  very 
similar  to  that  of  this  observation.  Therefore,  the  effect  of  the 
primary  beam’s  shape  on  the  radial  distribution  of  sources 
should  be  the  same  for  these  two  data  sets.  As  a  result,  the 
XMM-LSS  325  MHz  data  provide  a  good  background  radial 
distribution  for  this  comparison.  The  XMM-LSS  radio  survey 
also  observed  a  “blank”  field,  so  the  radial  distribution  of  the 
sources  is  dominated  by  instrumental  effects  and  not  by 
structure  intrinsic  to  the  program  source. 

We  also  compared  the  radial  distribution  of  sources  pro¬ 
jected  into  the  plane  defined  by  the  inclination  of  M31,  which 
may  contain  more  information  about  M3 1  itself  This  was  done 
by  first  converting  the  R.A.  and  decl.  of  each  source  in  the 
GLG  and  XMM-LSS  source  list  to  angular  coordinates  {X,  Y) 
relative  to  the  center  of  M3 1  (or,  for  the  XMM-LSS  sources, 
the  pointing  center  of  that  observation)  using  the  IDL  task 
WCSSPH2XY  and  then  translating  these  angular  coordinates 
into  the  plane  of  M31  on  the  sky  (Ymsi,  Tm3i)  using  the  fol¬ 
lowing  formulae: 

YM3i  =  Ycos52°-FTsin52°,  (1) 

Tm3i=  Tcos52°  -  Xsin52°,  (2) 

where  52°  is  the  angle  of  the  optical  disk  of  M31  relative  to 
north  on  the  sky  (Braun  1990).  Since  the  plane  of  M31  is 
inclined  77°  to  the  plane  of  the  sky,  a  circle  of  radius  R  in  M3 1 
would  be  an  ellipse  in  the  {X^^i,  Im3i)  coordinate  system. 
Using  this  fact,  the  distance  of  a  source  from  the  center  of  the 
observation  in  the  plane  of  M31  (Rm3i)  was  calculated  using 
the  following  formula: 


We  analyzed  the  statistical  properties  of  this  population  to 
see  if  radio  sources  in  the  field  of  M31  field  differed  from 


As  shown  in  Figure  6,  there  are  substantial  differences 
between  the  radial  distribution  of  the  XMM-LSS  and  GLG 
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Fig.  6. — Radial  distribution  {top)  and  i?M3i  distribution  {bottom',  defined  in 
§  4.1)  of  GLG  {filled  circles)  andXMAf-LSS  survey  {solid  line)  sources.  The 
diamonds  indicate  all  nonextended  GLG  sources,  while  the  triangles  are  non- 
extended  GLG  sources  with  5325  >  4  mJy  within  1?34  of  the  pointing  center,  i.e., 
the  limits  of  the  XMM-LSS  survey  (Cohen  et  al.  2003).  The  dot-dashed  line  is 
the  distribution  of  the  XMM-LSS  survey  scaled  by  the  ratio  of  “triangle”  GLG 
sources  to  XMM-LSS  sources  (see  §  4.1  for  details).  The  black  vertical  line  in 
the  bottom  graph  marks  the  outer  edge  of  the  optical  disk  of  M31. 

surveys.  The  first  ditFerence  is  the  number  of  sources  in  the 
XMM-LSS  (256)  and  GLG  survey  (405).  Some  of  this  differ¬ 
ence  is  due  to  differences  in  the  data  analysis — the  XMM-LSS 
list  does  not  include  sources  found  in  smoothed  maps  (“Ex” 
GLG  sources),  only  goes  out  to  l!34  from  the  pointing  center 
(in  the  GLG  survey,  the  furthest  source  is  1  ?75  from  the  pointing 
center),  and  does  not  include  sources  with  8^25  <  4  mJy  (Cohen 
et  al.  2003),  while  the  GLG  survey  has  a  5  ct  sensitivity  of 
^3  mJy.  However,  after  applying  the  above  limitations  to  the 
GLG  source  list  there  are  still  significantly  (2.9  a)  more  GLG 
sources  (302)  than  XMM-LSS  sources  (256).^  Some  of  this 
difference  may  be  due  to  the  fact  that  AMM-LSS  source  list  was 
derived  using  the  AlPS  task  VSAD  (Cohen  et  al.  2003),  a  task 
similar  to  SFIND  without  its  FDR  capabilities,  since  the  FDR 
algorithm  allows  SFIND  to  detect  fainter  sources  than  it  would 
without  it. 

^  This  implies  a  scaling  factor  of  302/256  =  1.18  between  the  two  sur- 
veys,  used  to  generate  the  dot-dashed  line  in  Fig.  6. 


The  second  difference  is  in  the  distribution  of  sources  be¬ 
tween  the  GLG  and  AMM-LSS  surveys — both  in  the  plane  of 
the  sky  and  in  Rm31-  In  fhe  plane  of  the  sky,  the  distribution  of 
the  AMM-LSS  sources  is  relatively  flat  between  0?25  and  1?25, 
while  the  distribution  of  the  GLG  sources  is  rising  in  this  range. 
This  difference  is  even  more  pronounced  in  the  Rm31  distri¬ 
bution,  where  for  Rm31~50  kpc,  there  are  noticeably  more 
GLG  sources  than  XMM-LSS  sources.  Inside  the  optical  disk 
of  M31  (Rm31  <  27  kpc),  there  are  145  GLG  sources  compared 
to  130  XMM-LSS  sources,  a  ~0.5  a  excess.  If  one  applies  the 
criteria  of  the  XMM-LSS  survey  to  the  GLG  source  list,  there 
are  in  fact  fewer  (119)  GLG  sources  than  AMM-LSS  sources. 
As  a  result,  we  expect  that  a  majority  of  the  GLG  sources  in  the 
optical  disk  of  M31  are  background  radio  sources. 

4.2.  Flux  Distribution 

If  the  325  MHz  population  of  M31  had  a  characteristic  flux 
density,  then  one  would  see  an  overdensity  of  GLG  sources  in 
a  particular  flux  density  range  when  compared  with  the  flux 
distribution  of  the  WENSS  and  XMM-LSS  surveys.  As  with 
the  radial  distribution,  the  primary  beam  affects  the  measured 
flux  distribution,  but  only  at  low  flux  densities.  Therefore,  we 
can  compare  the  flux  distribution  of  the  GLG  survey  to  that  of 
both  the  WENSS  and  AMM-LSS  surveys.  Figure  7  shows  that 
the  flux  distribution  of  GLG  sources  does  not  have  an  over¬ 
density  in  a  particular  flux  density  bin  or  range  of  flux  density 
bins — implying  that  the  flux  distribution  of  sources  in  M3 1  is 
not  substantially  different  from  that  of  background  sources. 
However,  given  that  we  only  expect  to  see  a  few  sources  in 
M3 1 ,  they  would  need  to  be  concentrated  in  one  or  two  flux 
density  bins  for  a  difference  fo  be  detected. 

In  order  to  determine  the  completeness  of  the  GLG  survey, 
we  compared  the  flux  distribution  with  a  model  determined 
from  the  327  MHz  WSRT  observations  (Wieringa  1991).  As 
shown  in  Figure  7,  there  is  good  agreement  between  the  model 
and  the  observed  flux  distribution  down  to  ^3  mJy — below  the 
5  a  sensitivity  of  the  observations  (4  mJy;  Wieringa  1991). 
Therefore,  we  estimate  that  the  GLG  survey  is  complete  to 
between  3  and  4  mJy. 

4.3.  Spectral  Properties 

In  general,  the  radio  spectrum  of  a  source  can  be  expressed 
as  a  power  law  attenuated  by  some  process.  In  this  paper,  we 
parameterize  the  attenuation  as  if  it  were  the  result  of  free- 
free  absorption  because  this  the  major  cause  of  low-frequency 
spectral  turnover  in  Galactic  objects.  Therefore,  we  fit  the 
observed  spectrum  to  the  following  function: 

where  Si,  is  the  flux  density  of  the  source  at  a  frequency  v,  S^o  is 
the  flux  density  of  the  source  at  frequency  uq  =  324.5  MHz,  a 
is  the  radio  spectral  index,  and  is  a  parameter  that  measures 
spectral  curvature.  If  the  spectral  curvature  is  truly  the  result  of 
free-free  absorption,  then  p  >  0  and  corresponds  to  the  optical 
depth  owing  to  free-free  absorption  at  =  1  GHz,  equal  to 
(Condon  1992) 

(^  =  3.3x10XEM)(^^)  ,  (5) 


EM  = 


(6) 
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Fig.  7. — Top:  The  flux  distribution  of  the  GLG  survey  {plus  signs),  WENSS 
survey  {solid  line),  and  the  XMM-LSS  survey  {dot-dashed  line).  The  y-axis  is 
SyMS),  where  n{S)  -  [7V(5o  <  5  <  5i)/(5i  -  5o)]/FOV  and  N{So  <  S  <  Sy) 
is  the  total  number  of  sources  with  flux  density  S  between  and  Si .  The 
term  corrects  for  Euclidean  geometry  (Condon  et  al.  1998).  The  diamonds 
represent  nonextended  sources.  The  vertical  lines  represent  the  completeness 
limit  of  theXMM-LSS  {left,  Cohen  et  at.  2003)  and  WENSS  {right,  Rengelink 
et  al.  1997)  surveys.  Bottom:  The  flux  distribution  of  the  GLG  survey  {plus 
signs)  compared  with  the  327  MHz  flux  distribution  model  {line)  derived  from 
WSRT  observations  (Wieringa  1991). 


where  EM  is  the  emission  measure  in  units  of  cm^®  pc,  Tg  is 
the  electron  temperature  (in  all  EM  calculations  presented 
here  we  assume  Tg  =  10"*  K),  rie  is  the  electron  density  in  units 
of  cm^^,  and  the  integral  in  equation  (6)  extends  over  the  line  of 
sight  to  the  object.  For  sources  with  an  insufficient  number  of 
counterparts  («ctrps  <  3)  at  other  radio  frequencies,  the  model 
used  does  not  include  the  attenuation  term  and  the  spectrum 
was  fit  to  a  simple  power  law. 

The  surveys  used  to  calculate  the  radio  spectral  index  (a) 
and  spectral  turnover  {(f)  of  the  GLG  sources  were  the  36 W 
{v  =  610  MHz;  Bystedt  et  al.  1984),  37W  {v  =  1412  MHz; 
Walterbos  et  al.  1985),  Braun  {v  =  1465  MHz;  Braun  1990), 
and  NVSS  {v  =  1400  MHz;  Condon  et  al.  1998)  surveys.  The 
properties  of  these  surveys  are  shown  in  Table  1 .  If  a  source  in 
another  catalog  was  less  than 


r  =  Max 


(7) 


away  from  a  GLG  source,  they  were  classified  as  a  match, 
where  is  the  major  axis  of  the  GLG  source,  is  the 
minor  axis  of  the  GLG  source  (these  are  not  the  deconvolved 
values  but  the  size  of  the  source  in  our  images),  9m  is  the 
deconvolved  major  axis  of  the  catalog  source,  and  9^  is  the 
deconvolved  minor  axis  of  the  catalog  source  (Braun  1990).® 
A  “good  match”  between  a  GLG  source  and  a  source  in  one  of 
these  catalogs  was  when  the  GLG  had  only  one  match  in  the 
other  catalog  and  its  counterpart  only  had  no  other  matches  in 
the  GLG  source  list.  The  results  of  this  comparison  are  sum¬ 
marized  in  Table  4.  The  number  of  false  matches  was  calcu¬ 
lated  by  shifting  the  position  of  the  radio  catalog  in  question  by 
every  combination  of  ±1',  0'  in  R.A.  and  ±1',  0'  in  decl. 
(except  for  a  shift  of  [0,  0]),  and  calculating  the  number  of 
matches  and  good  matches  in  each  combination.  The  number 
of  false  (good)  counterparts  reported  in  Table  4  is  the  average 
number  of  (good)  matches  detected  in  the  eight  shifted  data 
sets,  and  the  error  in  these  values  is  the  standard  deviation  of 
the  number  of  (good)  matches.  As  seen  in  Table  4,  the  number 
of  false  counterparts  and  false  good  counterparts  are  much  less 
that  observed.  The  WENSS  survey  has  a  larger  number  of  false 
comparisons  than  the  other  surveys  because  its  beam  is  much 
larger,  resulting  in  a  larger  search  radius.  In  order  to  avoid 
problems  with  source  confusion,  a  and  (p  were  calculated 
using  only  values  from  a  “good  match”  by  fitting  the  spectrum 
to  the  form  defined  in  equation  (4)  using  a  least  mean  squares 
algorithm. 

For  GLG  sources  without  a  good  counterpart  in  any  of  these 
radio  surveys,  we  calculated  the  upper  limit  on  the  spectral 
index  a .  In  order  to  do  this,  we  made  a  sensitivity  map^  of 
the  36W,  37W,  and  Braun  surveys  (shown  in  Fig.  1)  by  first 
making  a  sensitivity  map  of  each  pointing  in  the  survey  and 
then  combining  them  together,  weighing  each  pointing  by  l/tr^, 
where  cr  is  the  rms  of  the  pointing  map.  The  upper  limit  was 
then  determined  by  adjusting  the  stated  minimum  detectable 
flux  density  by  the  ratio  of  the  survey’s  sensitivity  at  the  lo¬ 
cation  of  the  source  to  the  survey’s  peak  sensitivity.  This  was 
not  necessary  for  the  NVSS  survey  because  our  FOV  is  entirely 
within  theirs,  and  an  upper  limit  of  5i4oo  =15  mJy  (the  com¬ 
pleteness  limit  of  the  NVSS  survey)  was  used  (Condon  et  al. 
1998).  The  spectral  index  of  these  sources  given  in  Table  3 
is  the  flattest  spectrum  consistent  with  all  nondetections.  The 
spectral  energy  diagram  (SEDs)  of  all  GLG  sources  can  be 
found  in  Figure  8,  and  the  calculated  a  and  f  of  each  source 
can  be  found  in  Table  3. 

For  sources  with  a  “good  match”  in  the  36W  (610  MHz) 
survey  and  in  at  least  one  of  the  1400  MHz  surveys  (NVSS, 
37W,  or  Braun),  it  was  possible  to  compare  their  a  32®  and 
a 325°  spectral  indices,  the  spectral  index  between  325  and 
610  MHz  and  325  and  1400  MHz,  respectively.  For  sources 
with  counterparts  in  multiple  1400  MHz  radio  catalogs,  a  32°° 
is  the  weighted  average  of  the  individual  spectral  indices. 
Figure  9  shows  the  graph  of  a  32°°  versus  afjs  for  the  156  GLG 
sources  that  met  this  criteria.  These  sources  are  not  evenly 
distributed  around  the  0325°  =  <3325  line;  99  GLG  sources  have 
a  325  >  0:325°,  while  only  57  have  032°  <  032°°.  A  larger 
number  of  032°  >  0^2°°  sources  is  expected  because  the  syn¬ 
thesized  beam  of  the  36 W  survey  is  larger  than  that  of  any  of 


^  Because  the  36W  and  37W  surveys  were  done  with  the  Westerbork  Radio 
Telescope  (WSRT)  as  opposed  to  the  VLA,  there  is  a  systematic  offset  in 
position  between  them  and  the  GLG  survey.  However,  this  difference  is  , 
much  smaller  than  the  >6"  radius  used  to  search  for  counterparts. 

^  A  sensitivity  map  shows  the  rms  of  an  observation  over  the  field  of  view. 
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TABLE  4 

Comparison  with  Other  Radio  Surveys 

Number  of  GLG  Sources  Number  of  GLG  Sources  Number  of  GLG  Sources  Number  of  False  Number  of  False 


Survey  in  Survey  FOV  with  Counterparts  with  Good  Counterparts  Counterparts  Good  Counterparts  Reference 


WENSS....  405  121  89  35.4  ±  8.2  16.3  ±  4.5  Rengelink  et  al.  (1997) 

36W .  357  176  161  8.8  ±  2.1  4.3  d=  1.8  Bystedt  et  al.  (1984) 

37W .  227  153  145  1.8  ±  1.6  1.4  ±  1.3  Walterbos  et  al.  (1985) 

Braun .  87  67  53  9.8  ±  2.6  4.9  ±  2.4  Braun  (1990) 

NVSS .  405  189  185  2.3  ±  1.6  2.1  ±  1.6  Condon  et  al.  (1998) 


the  1400  MHz  surveys,  meaning  that  the  reported  flux  density 
of  36W  sources  possibly  include  emission  resolved  out  in  the 
1400  MHz  surveys.  However,  23  of  these  sources  have  |  0325“ 
0325®  |>  3Max((Ta6i»,  cr„i4oo)^  where  is  the  error  of  ajjj 
and  CTaWoo  is  the  error  of  a  statistically  significant  dif¬ 

ference  that  cannot  be  explained  by  a  difference  in  beam  size. 
All  but  one  of  these  sources  have  ctjj®  >  oi\^^ — the  signature 
of  both  free-free  absorption  and  synchrotron  self-absorption, 
and  such  sources  have  been  observed  in  previous  radio  surveys, 
e.g.,  De  Breuck  et  al.  (2000).  GLG  186,  the  only  such  source 
with  01325  <  01325° — has  been  categorized  a  variable  source 
(see  J.  D.  Gelfand,  T.  J.  W.  Lazio,  &  B.  M.  Gaensler  2004,  in 
preparation,  hereafter  Paper  II,  for  more  information). 

To  determine  if  the  GLG  sources — both  in  and  out  of 
the  optical  disk  of  M31 — had  different  spectral  properties 
than  background  sources,  we  compared  the  ^325-0  distribution 
(Zhang  et  al.  2003)  and  the  o  number  distribution  (Figs.  7  and 
8  in  De  Breuck  et  al.  2000)  of  the  GLG  sources — shown  in 
Figure  10 — with  that  of  sources  in  both  the  WENSS  and  NVSS 
survey.  Both  the  5325-0;  and  a  number  distribution  of  GLG 
sources  is  consistent  with  the  WENSS/NVSS  results.  We  also 
compared  the  spectral  index  distribution  of  GLG  sources  with  a 
higher  frequency  radio  counterpart  inside  and  outside  the  op¬ 
tical  disk  of  M3 1  to  see  if  there  are  any  differences  between 
these  populations,  also  shown  in  Figure  10.  While  there  are 
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Fig.  8. — Radio  spectrum  (u  —  325—1400  MHz)  of  all  GLG  sources.  A  filled 
circle  represents  a  good  match,  an  arrow  represents  an  upper  limit,  an  open 
square  on  the  jc-axis  represents  a  “bad  match”  (the  GLG  source  either  had  more 
than  one  counterpart  in  this  catalog  or  the  catalog  source  had  more  than  one 
counterpart  in  the  GLG  source  list),  and  the  dashed  line  is  the  best  fit  to  eq.  (4). 
The  upper  limits  are  scaled  by  the  local  rms  of  the  survey  in  question.  [See  the 
electronic  edition  of  the  Journal  for  additional  panels  of  this  figure.] 
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substantial  selection  effects  involved,  e.g.,  the  FOV  of  the 
deeper  1 .4  GHz  surveys  are  concentrated  on  the  optical  disk  of 
M3 1  so  steep  sources  outside  M3 1  might  be  missed,  there  is  no 
noticeable  difference  in  the  number  distribution  of  a  between 
these  regions.  However,  this  is  not  unexpected  since  the 
dominant  radio  sources  of  a  galaxy — SNRs  and  H  ii  regions — 
have  similar  spectral  indices  distributions  than  that  of  the 
dominant  forms  of  background  sources — FR  I  and  FR  II  radio 
galaxies. 


5.  CONCLUSIONS 

In  this  paper,  we  present  the  results  of  a  high-resolution  (6") 
325  MHz  survey  of  a  7.6  deg^  field  centered  on  M31.  In  this 
data  set,  we  identified  405  discrete  sources  (the  “GLG”  source 
list)  and  determined  their  size,  orientation,  and  flux  density. 
Through  comparisons  with  610  and  1400  MHz  radio  surveys 
of  the  same  region,  we  calculated  the  spectral  index  a  and 
spectral  curvature  (f  of  these  sources. 

To  determine  the  nature  of  the  GLG  source  list,  we  com¬ 
pared  its  properties  to  that  of  source  lists  derived  from  blank- 
field  and/or  large-area  radio  surveys — the  WENSS,  NVSS,  and 
XMM-LSS  radio  surveys.  A  comparison  of  the  radial  distri¬ 
bution  of  GLG  and  XMM-LSS  radio  sources — both  in  the 
plane  of  the  sky  and  in  the  plane  of  M3 1 — showed  that  there  is 
not  a  significant  concentration  of  sources  in  the  optical  disk  of 
M31.  However,  there  does  appear  to  be  an  excess  of  GLG 
sources  at  the  edge  of  the  field  of  view.  The  reason  for  this  is 
unknown,  but  most  likely  caused  by  the  different  source 
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Fig.  9. — “Color-color”  diagram  of  all  156  GLG  sources  with  a  good  match 
in  the  36W  catalog  and  at  least  one  of  the  1400  MHz  catalogs.  The  diamond  is 


the  one  extended  source  that  met  this  criteria,  and  the  squares  are  sources  with 
«32S  -  “32“  l>  3Max((7„6,o,o-„,4»). 
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Fig.  10. — Clockwise  from  top  left:  a  vs.  5*325  for  the  244  GLG  sources  with  a  good  counterpart  in  at  least  one  other  radio  catalog;  median  spectral  index  vs. 
325  MHz  flux  density  for  these  sources  (data  points  are  a  binned  version  of  the  top  left  graph,  and  the  error  bars  are  the  standard  deviation  of  spectral  index  in  the 
5*325  bin)  overlaid  with  the  results  from  the  WENSS-NVSS  surveys  {triangles',  Zhang  et  al.  2003);  fraction  distribution  of  a  for  sources  inside  {solid  line)  and  outside 
{dashed  line)  the  optical  disk  of  M31;  and  number  distribution  of  a. 


finding  programs  used  between  the  surveys.  A  comparison 
between  the  flux  distribution  of  GLG  sources  with  that  of 
the  WENSS  and  AMM-LSS  radio  surveys  revealed  no  statis¬ 
tically  signiflcant  overdensity  or  underdensities  in  a  sizable 
flux  range.  The  a  number  distribution  of  GLG  sources  was 
similar  to  that  of  sources  in  both  the  WENSS  and  NVSS  radio 
surveys,  as  was  the  relationship  between  325  MHz  flux  density 
5325  and  a.  The  spectral  index  distribution  of  sources  within 
the  optical  disk  of  M3 1  was  also  similar  to  that  of  sources  in 
outside  the  optical  disk  of  M3 1 .  All  of  this  implies  that  the 
GLG  source  list  is  dominated  by  background  radio  galaxies 
and  not  sources  intrinsic  to  M3 1 . 

In  an  accompanying  paper  (J.  D.  Gelfand,  T.  J.  W.  Lazio,  & 
B.  M.  Gaensler  2004,  in  preparation),  we  compare  the  GLG 
source  list  with  that  of  far-IR,  IR,  optical,  and  X-ray  catalogs 
of  this  region  and  use  this  multiwavelength  information,  as 
well  as  the  morphology  and  spectral  index /curvature  of  a 
source,  to  determine  the  nature  of  individual  sources.  Through 
the  method  described  in  that  paper,  we  have  identified  five 
supernova  remnant  and  three  pulsar  wind  nebulae  candidates 
in  M3 1 ,  as  well  as  a  galaxy  merger,  BL  Lac  candidate.  Giant 


Radio  Galaxy  candidate,  and  several  low-frequency  variable 
and  ultrasteep  spectrum  (a  <  — 1.6)  sources. 
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